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A model is presented in which a single scalar field is responsible for both primordial inflation 
at early times and then dark energy at late times. This field is coupled to a second scalar field 
which becomes unstable and starts to oscillate after primordial inflation, thus driving a reheating 
phase that can create a high post-inflation temperature. This model easily avoids overproduction of 
gravity waves, which is a problem in the original quintessential inflation model in which reheating 
occurs via gravitational particle production. 
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I. INTRODUCTION 

Peebles and Vilenkin proposed in [l[ that both inflation 
and dark energy could be a result of the same scalar held, 
with vacuum expectation value </>, interacting only with 
gravity and itself via the potential term V{4 1 ), which they 
chose to be 



V (</>) = A(0 4 + Af 4 ) 
AM 8 
6 i + M A 



for <f) < 
for cb > 0. 
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At tree level the evolution of the vacuum expectation 
value of the </>-field is governed by 



dV 



(2) 



with the cosmological expansion rate related to energy 
density p by the Friedmann equation, 



= [ - 
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where the Planck mass m p i = G" 1 / 2 = 1.22 x 10 19 GeV 
and overdots signify derivatives with respect to coordi- 
nate time. 

In the original scenario the universe begins domi- 
nated by the potential energy of the scalar field 4>. The 
field has some large, negative value, and slowly rolls to- 
ward the origin. The slow-roll conditions of inflation are 
satisfied, and the universe expands exponentially. During 
inflation, fluctuations in 4> are frozen into the field. These 
seed future structure in the universe. To provide the cor- 
rect level of fluctuations, as usual A ~ 10~ 14 . When 
4> ~ — m p i, the inflationary epoch draws to a close. Thus 
far the situation is identical to <p 4 chaotic inflation Q, 
but from here on it differs. The kinetic energy ^<fi 2 of the 
field is no longer negligible, and it soon begins to dom- 
inate the universe [3j . This phase is termed "kination" . 
The field behaves approximately as stiff matter, with its 



energy density redshifting as cx a -6 . Taking the po- 
tential energy at the end of inflation to be the same as 
the kinetic energy, a simple estimate gives 
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with solution, a cx i 1 / 3 and thus t = (3-ff) -1 . Here the 
subscript x indicates the value at the end of inflation. 
During this kinetic dominated phase, the field moves as 



-m p iln (a/a x ) - m pl . 



(5) 



Thus, spacetime started in approximately de Sitter form, 
and ended dominated by the kinetic energy of a homo- 
geneous field. Ford showed [3] this transition leads to 
gravitational particle production. The result is a small 
but important energy density of relativistic particles, 



(G) 



where N s is the number of scalar fields. Thcrmalization 
was found in [l[ to occur at a radiation temperature 



T th ^ 10 9 iV 3/4 GeV . 



(7) 



The radiation redshifts away slower than the energy in 
the field. Provided M in Eq. ^ is not too large, Pee- 
bles and Vilenkin showed the universe transitions to a 
radiation dominated epoch, with a temperature of 



r, 



R.H 



10 3 7V 3 / 4 GeV . 
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The field remains essentially static from this point on, 
mimicking a cosmological constant. The value of M can 
be then chosen so that V(4> r ) matches today's observed 
value of dark energy density. In the analysis in this 
turned out to be M - 10 6 GeV. 

Unlike most inflationary scenarios, in the one by Pee- 
bles and Vilenkin [1] at the end of inflation the inflaton 
field does not undergo a series of damped oscillations 
about a potential minima, which is the mechanism by 
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which reheating usually takes place Q. The absence of 
this behaviour means gravitational particle production 
has to be relied upon instead. The same gravitational 
mechanism works to produce a stochastic background of 
gravitational waves (GW) [|0,1, and the overproduc- 
tion of gravity waves is one of the potential dangers of this 
model. Gravitons behave as minimally massless scalar 
fields, and the energy density of the gravitons by the end 
of inflation is just that of a single scalar field, times two 
polarization states, so that the ratio of the energy densi- 
ties in GW to radiation at the end of inflation is simply 
{PGw/Pt)x — *2/N s . At the time of nucleosynthesis paw 
will contribute as an effective extra degree of freedom 
in radiation, but the success of Big Bang Nucleosynthe- 
sis (BBN) in predicting the abundances of light elements 
puts a constraint on the extra number of light degrees 
of freedom, be those neutrinos species or gravitons, such 
that (paw I Pt)bbn < 0.2; tracing back the evolution of 
this ratio upto the start of the kination period this would 
sets a lower bound N s > 100 [J. 

The problem of gravity waves overproduction is quite 
generic in models of inflation followed by a long period 
of kination (9[ , or in general stiff matter domination [l(| , 
like brane world inflation [ll[ . To allow for a more effec- 
tive reheating process, able to suppress the relative con- 
tribution of the GW at the time of BBN, one can invoke 
alternative methods like instant preheati ng Il2l| , curvaton 
reheating '13|, or Born-Infeld reheating |l4j. Any alter- 
native implies introducing extra scalar degrees of freedom 
at the time of inflation, like in curvaton reheating, and/or 
direct couplings of the inflaton field to the light degrees 
of freedom as in instant preheating. 

In this letter we explore the possibility of a simpler 
scenario, recovering a more typical reheating mechanism 
driven by the decay of an oscillating massive field Q. 
We extend Peebles and Vilcnkin model by a new scalar 
field x coupled to the inflaton field [l5[ with a hybrid- 
like potential [l6| . In our scenario, once the inflaton field 
falls below a critical value, the x ne ld can start oscillat- 
ing, thus gaining energy that afterwards can be converted 
into radiation through pcrturbative decay as in the usual 
reheating mechanism. However, in contrast to the stan- 
dard reheating picture, in our scenario reheating takes 
place during kination instead of the more standard mat- 
ter domination. Unless the perturbative decay of the x 
field is tiny, this results in a larger reheating T than in the 
original model of Peebles and Vilenkin [l| , and a shorter 
kination phase. The more efficient reheating also ensures 
that radiation domination takes over kination, well be- 
fore the inflaton vacuum energy starts dominating again. 

This letter is organized as follows. In Section II the 
potential and parameters of the model are set. Also the 
general behaviour is described of the field x after infla- 
tion, when it can oscillate and drive reheating. The re- 
heating temperature Trh is then computed in Section 
III. For the mechanism to work, we need to check first 
that x indeed oscillates when the inflaton field passes 
through the critical point, and that it does not backreact 



on the evolution of 4>. Fulfilling these conditions sets the 
constraints on the model parameters, which are given in 
Section IV. In Section V we present the range of Trh 
consistent with the constraints. The transition to a radi- 
ation dominated universe is studied in Section VI. Once 
the constraints are fulfilled, the transition before the on- 
set of dark energy domination is practically ensured. As 
this is a hybrid-like model, we comment on the issue of 
domain walls in section VII. Finally in Section VIII we 
present the summary and future work related to quantum 
corrections. 



II. GENERAL BEHAVIOUR OF THE HYBRID 
FIELD 

To enable once more the standard reheating, we intro- 
duce a new scalar field x, that we couple to the inflaton 
field. The effective potential at tree-level is 

U(cf>,x) = V(cb) + ^x 2 ^ 2 ~rn 2 ) + ^x\ (9) 

where V((j>) is the potential given by equation fl}. The 
parameters g, to, and X x are as yet undetermined con- 
stants. Note that when x is relaxed near the origin, 
U(<f>, X )*V(<t>). 

The x-field is assumed to be located somewhere near 
its minima at x — 0. As the <f> field evolves from large 
negative values to large positive ones, the turning point 
at the origin temporarily becomes unstable for the x~ 
field, and two minima are generated on either side. The 
position of the minima, Xmin are given by 

Xlsn = f K - 2 ) . (10) 

They exist only while \(f>\ < to. The bottom of the well is 
at a negative value of potential energy. This is a result of 
how the energy of the system has been defined. It has no 
physical consequence for the present scenario, since the 
total energy density will always remain positive, as it is 
dominated by the kinetic inflaton energy density. 

The x-field has a characteristic response time r to react 
to changes in the potential. This can be estimated as 
t _1 ~ rn x {4>) ~ ^(to 2 — 2 ) 1 / 2 , provided <j> is not too close 
to m. If changes in the potential take place quicker than 
this response time, the field will have no chance to react. 
For instance, if <f> moves between — m and +to in a time 
At <C t throughout, x w iU have no time to move before 
the origin becomes a stable minima once more. On the 
other hand, if the changes take place on timescales that 
are much longer than the response time, i.e. At 3> t, the 
X field will be able to relax into the minima very quickly. 
If this is the case throughout, there will hardly be any 
oscillations (and hardly any reheating). 

A fact alleviates the above difficulties: (j> is slowing 
down, and is doing so quickly (its evolution is logarithmic 
in time). This means the timescale At becomes gradually 
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longer. Furthermore, the timescale r that governs how 
quickly \ reacts is not static. It is longest when \<f>\ w m 
and is shortest at cj) = 0. We will take the view that the 
field begins to move before <fi sa 0. This places an imme- 
diate constraint on the model parameters: At 3> t|0=o- 
We should immediately note the better this inequality is 
satisfied, the sooner x will start to move after becoming 
unstable. This will limit the amount of energy for its 
oscillations. 

As a simplified picture, consider the x field to be frozen 
at the origin up until a time r after becoming unstable. 
Afterwards, the potential will be approximately constant 
during the fast oscillations of the field. The energy in 
the oscillations (available for reheating) will then be the 
depth of the potential well at the point when the field 
begins to move. This will be some fraction f(g, m) of the 
maximum potential well depth. The energy available for 
reheating can then be written: 



JO) 



f(g,m) 



4 4 

m g 
4A V 
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Eventually, after many oscillations, <f> reaches +m, and 
the x ~ field then oscillates about its stable minima at 
the origin. These oscillations will reheat the universe. 

The fraction f(g,m) now needs to be estimated. For 
this, we assume the x -field moves some short time after 
it becomes unstable, with the instability occurring at <f> = 
—m. As this time signals the start of reheating, it will be 
called t xe . Writing the time since the start of reheating 
as St, then we can estimate that the x field will begin 
to move when r ~ St. In other words, the field moves 
after it has been unstable for a time approximately the 
same as its characteristic reaction time (which is itself a 
function of St). The depth of the well is given by 



4 4 

TT -4/2 /2\ 2 m 9 
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Comparison with equation shows 

f(g, m ) = m ~ 4 i™ 2 - (f 2 ) 



(12) 



(13) 



with <f> evaluated when St ~ t. We assume <p moves only 
slightly, an amount Scj) — 4> — m, before x begins to move. 
Re-writing equation (f5]) gives 



(14) 



We wish to find the value of St satisfying St ~ r. As 
t _1 ~ g(m 2 — cj) 2 ) 1 / 2 , we can insert cj) — Scj) ~ m to find 



-2 2 2 / 6 St 

T ~ m9 3 V 

The same substitution into f(g,m) yields 

WO 2 



f(g,m) 



(15) 



(16) 



to leading order in Scj). 

Applying the condition St ~ r, we can find the value 
of St when the x ~ field begins to move. Relating Scj) to 
St with equation (|14[) . we finally obtain 

/( 5 ,-)-(^) 2/3 (^f) 8/3 (-pAe)- 4/3 5 - /3 . d7) 

Note that this order of magnitude approximation breaks 
down if x begins moving when 5(f) is not small compared 
to m, signalled by /(g,m) approaching (or exceeding) 
unity. This expression should not be trusted in such 
circumstances. Assuming Scj) to be small is a slightly 
stronger constraint than the constraint discussed earlier, 
that At > r. The former gives the condition the field 
moves quickly after becoming unstable, the latter is just 
the condition the field moves at all, at or before ^0. 
We will return to this point when discussing the param- 
eter constraints in Section IV. 



III. REHEATING TEMPERATURE 

Once cj> > m, the x field will return to the origin. If 
it acquired kinetic energy due to its temporary displace- 
ment, it will now oscillate about the origin. The previous 
section established the amount of energy expected from 
these oscillations. 

If the oscillations are small, the field can be approxi- 
mated as undergoing simple harmonic motion. In such a 
case, the elementary theory of reheating can be applied 
[5(. In this phenomenological approach, an extra term 
T x x is added to the equation of motion of the field to 
account for particle decay. The value of T x is taken to 
be the decay rate of the particle. The field x then obeys 
the equation of motion 



x- 



3-X 

a 



dU 

dx 



(18) 



This approach is valid when the oscillations are small and 
the oscillations are well approximated as simple harmonic 
motion. For a more complete picture, valid at the early 
stages we should also consider the effect of preheating [17| 
(see also [l2j for an analysis of preheating in the context 
of quintessential inflation). However, these details will 
be ignored in this paper and we will examine only the 
simplest reheating estimates. 

We stress that this phenomenological approach is only 
valid while the field is undergoing coherent oscillations 
about its minima. It will not be valid otherwise, nor 
over short timescales, and is likely to fail when \(j>\ ~ 
m. For this reason we will consider reheating to take 
place only while <p > m, so that we may have confidence 
our calculations are always carried out in an appropriate 
regime. 

We assume H, <fi and c\> can be taken to be approxi- 
mately constant over a single oscillation. Equation [)18p 
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is re-written by replacing x 2 by its value over a com- 
plete oscillation, (x 2 ) cycle — Pxi which is valid for simple 
harmonic motion. This yields 



Px + % H Px + T xPx = 9 2 {X 2 )H ■ 



(19) 



Assuming the field is still undergoing simple harmonic 
\p x = V(cp, x) = \g 2( t )2 {x 2 )- Then we can write 



motion, ^p x 



Px + ' iH Px + T xPx = tPx ■ 



This can be solved, 



Px 



P ( x l) 



V a 



0(t) 



exp [-r x (* - t m )] 



(20) 
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In 



(22) 



We have used subscript m to indicate the value of a vari- 
able when (f> — m. From energy conservation, it follows 
that the radiation density must obey 



p r + 4Hp T = Y xPx 



(23) 



Imposing the condition there is no radiation at the start 
of decay, an approximate solution is found by neglecting 
the exponential decay of p x . This will be valid up to 
t « r^ 1 . After this, the energy in the x field will decay 
rapidly away and the radiation will simply redshift with 



its usual 



behaviour. 



Inserting our earlier expression for p x into equation 
([23]) . the solution for the radiation energy density can be 
written as 



^X L X L m 



(l-(t/t m )- 4/3 )(l-b) 

+ &(4/3)ln(t/t m ) 



(24) 



with b = \\ ■§- and the constant of integration cho- 
sen so that there is no radiation at cf> = m. 

Assuming (i/t ro )~ 4 / 3 <C 1 before t reaches T" 1 , the 
energy density in radiation should be well approximated 
by 



-o {m) T t 



1-6+ -b In (t/t m ) 



(25) 



The energy density of the radiation continues growing 
logarithmically, despite the energy loss from redshifting. 
This is due to the mild amount of energy being added to 
the x field by its coupling to <f>. It will continue to grow 
in this way until t w T^ 1 . If this is a sufficiently late time 
that the last term in Eq. ([25]) dominates, then 



-o (m) r / 
p x 



Ht/t m ) . 



(26) 



FIG. 1: Numerical solution for the motion of the \ field, with 
j = 5x 10~ 4 , X x = 1, m = 5 X 10 16 GeV and a = 1 x 1(T 4 . 
The field begins at rest, with a small displacement from the 
origin. The resulting oscillations continue when the field re- 
turns to oscillate about the origin, and are well approximated 
as damped, simple harmonic motion. 



Once the universe has had chance to thermalize, the tem- 
perature is related to the energy density via, 
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(27) 



where g» is the number of degrees of freedom. 

Thermalization occurs when the interaction rate n r a 
becomes comparable to the expansion rate H, where n r is 
the number density of the light degrees of freedom and a 
the their interaction cross section. We can estimate that 
the light degrees of freedom are created with a typical en- 
ergy ui r~j pl^ 4 (a rc /a), and a ~ a g /uj 2 , with a g being the 
strength of the mediating interactions. Using Eq. (|2l)|) 
with the expression of t le ~ t m given in the next section 
Eq. (|31[) . and as a typical value for the coupling in the 
cross-section a g ~ 0.01, one gets that at the beginning 
of the reheating period, 



1/4 
Pr a g 



lO^lOa 1 / 4 ) ( — 

. TOpl 



1/3 /_ll/3\ 



(28) 

We have also written the decay rate for a massive particle 
as r x ~ am x , where a is the coupling constant mediat- 
ing the decay, and m x ~ gm the x mass. For the analyses 
of the reheating done in this section, the decay rate and 
therefore the coupling a must be such that T x t m <C 1. 
For values of the parameters consistent with the con- 
straints given in the following section, this is general the 
case with a ~ 10 -4 . Significantly, with this choice for a 
then Eq. (j2"5|) is already larger than one. Therefore, light 
degrees of freedom thermalized promptly after they are 
produced. 

We have confirmed these approximations are successful 
in their appropriate regimes by numerically solving this 
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system of differential equations (Friedmann's equation, 
the equations of motion for both fields, and the radiation 
energy density) . Figure [1] shows the numerically deter- 
mined evolution of the x ~ field as it becomes unstable, 
and the subsequent damped oscillations that reheat the 
universe. 



IV. PARAMETER CONSTRAINTS 

We have made two assumptions that can be formu- 
lated as simple constraints on combinations of parame- 
ters. The first of these is that \ begins to move toward its 
new equilibrium at or before <fi ~ 0. Earlier we noted this 
condition was At ^> r. The second assumption is that 
the x field does not significantly influence the motion of 
the (j> field. It could do this either from the coupling, or 
from the energy density of the field modifying the expan- 
sion rate of the universe. 

First we shall calculate At. As cf> — —m at t rc , we 
can use equation (|22l with these values inserted. Re- 
arranging, 



/487T- 



t —te 
and by writing At = t m — t Te , gives 



At = tre 



1 



(29) 
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with the rightmost expression in the case m/m p \ <C 1. 
Thus, it is noteworthy that t ro w t m is a good approxi- 
mation. 

An estimate of t re is still needed. We use equation ((3]) 



and H 
factors, 



\t 1 , to write it in terms of a ratio of scale- 



24?rAr 



(31) 



'pi 



where a ro is the scale factor at t 
move between —m and +m is then given by 

Atc(^) 3 



The time for <fr to 



2 m 
A ml 



(32) 



The ratio of scale- factors can be found from equation ([5]). 
Re-arranging and setting 4> — —m at a — a Ie , 

a ie /a x = e^- m ^V^ ~8. (33) 

As such the time interval available can be written simply 
as 



At ~ 1(T 9 (GeV)~ 

m p i 



(34) 



with A = 10 14 . Recalling At 3> r, a constraint on the 
model parameters can be constructed: 



rn 
m pl 



.9» 10 



-10 



(35) 



Now that we have an expression for t re , a constraint on 
f(g,m) -C 1 can also be calculated. Replacing Eq. (pITj) 
into (fTTf gives 



\m p ij \a rc J 



(36) 



This is comparable although slightly more restrictive 
than Eq. (j3"5|) . This suggests that there can be a region in 
parameter space where the second constraint is violated 
(suggesting 1 > f(g,m) > 0.1) but the first satisfied (so 
that the field still begins to move out of its unstable posi- 
tion) . The resulting reheating temperature in this region 
is also fairly insensitive to changes in g and m, compared 
with when the second constraint is well satisfied (and 
the reheat temperature influenced by f(g, m) as given by 
equation (fT7]0. Nevertheless, to ensure that we are in 
the region of parameter space for which the field has os- 
cillated enough to reheat the universe, when referring to 
these constraints we will use (m/m p \) 2 g > 10 -8 . Having 
m p i as the largest possible mass scale in the model , the 
strength of the <p~x interaction is therefore bounded from 
below with g > 10~ 8 . 

A second constraint exists from the requirement that 
the x field does not influence the motion of the (j> field. 
First we will consider the requirement the negative po- 
tential energy from the coupling term is not significant 
compared to the kinetic energy of 4>. As the kinetic en- 
ergy is constantly diminishing, it will be simpler to over- 
estimate the potential energy and underestimate the ki- 
netic energy. We will therefore take the potential energy 
to be its maximum value, and the kinetic energy to be 
the value at <f> = m. No matter the evolution of these 
quantities, if the quantities evaluated at these two time 
intervals are not comparable, they never will be. Read- 
ing off the kinetic energy from equation (|U) and requiring 
this always greatly exceeds the maximum of the negative 
potential energy gives, 



Am pl (^M 6 > m494 



or 
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If we choose values to satisfy At ^ r, 

^- < 10~ 3 . 
A v 



(37) 



(38) 



(39) 



Now consider the effect of the coupling on the motion 
of <fi directly. The equation of motion gives 



3H(f> 



2 2 i 

9 x <i> 



0. 



(40) 



where the kinetic energy is treated as the dominant con- 
tribution to energy of the 4> field. The third term is the 
contribution due to the coupling, and we wish to ensure 
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FIG. 2: The range for the reheating temperature Trh in the plane g — m/m p i, for X x = 1 (left-hand panel) and A x = 0.1 
(right-hand panel). Decreasing the value of A x has the effect of reducing the allowed region of parameter space. The shading 
density indicates high (light) to low (dark) reheating temperature in the allowed region of parameter space. 



this is negligible compared to the second. Using equation 
(|3"|) to express Hubble's parameter in terms of the kinetic 
energy ss l/2</> 2 , this condition can be written as, 



I An 6 ,2 
3 m p i 



(41) 



Treating the right-hand side as the average value over 
oscillations of the \ held, we can rewrite it using p x — 
g 2 4> 2 (x 2 )- Writing these energy densities out explicitly, 
including their evolution with scale-factor gives 



jM _£_ P -r x (t-t m ) 



— — »f(9,m) 



4A, 



(42) 



The right-hand side will quickly become tiny when t ~ 
r" 1 . We need only consider the value of a oc t 1 / 3 when 
(a m /a) 3 ~ T x t m . Using T x ~ agm 



this occurs, and then 
and Eq. (JTTJ) gives: 



(^) 2/3 ^«V2, 
V m ) Ay 



aA "l/6 ( ^ 



3xl0 3 a, (43) 



and thus unless the decay rate is tiny and the reheating 
period too long, once we fulfill the other constraints this 
is practically always fulfilled. In other words, having the 
parameter values such that the field x performs some os- 
cillations around the minimum while the expansion rate 
is dominated by the inflaton kinetic energy, those oscil- 
lations will not backreact onto the motion of the inflaton 
field. 



V. RANGE OF TEMPERATURES 

The constraints show there can be significant varia- 
tion in the resulting temperature. The energy density 
depends most sensitively on g, and the range this pa- 
rameter can take is severely constrained by the other 
parameters. If A x is too small, g must be made small 
enough to avoid interfering with the evolution of (p. If 
m is made too small, g must be made large enough to 
ensure the field reacts while it is unstable. For instance, 
if A x ~ l then g can range from g ~ 1CU 2 at its largest 
(with m/m pl - 10~ 3 ) to g ~ 1(T 6 (with m/m p i - KT 1 ), 
assuming m is kept below the Planck scale. Decreasing 
A x can constrain it further. 

T~ l . Plugging Eqs. 
agm and A = 10 - , we have: 



Reheating ends by the time t 
(EU, (HZD, with r 

4/3 

p r c± S x 11) "m^a 



(—) 

\m p ij 



11/3 



ln(T x t 



(44) 



1 and a 



Therefore, working for example with A x i 
10~ 4 , we find the range 

9 X J V RH ~ (10 11 - 10 14 ) GeV. 
One can also look at very weak coupling, where for exam- 



(45) 



pie g ~ 1(T 4 , A x 



io- 



10 , and m/m p i ~ 10 
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gives g* /4 T RH ~ 3 x 10 13 GeV. In Fig. we give the 
range for the reheating temperature Trh m the plane 
g — mj TOpi for two values of A x . As we decrease the value 
of the x self-coupling, the allowed region in the plane will 
be further reduced, and that would be the main effect on 

VI. TRANSITION TO RADIATION 
DOMINATION 

The kinetic energy of the <j) field is dropping quickly. 
Peebles and Vilenkin noted that unless radiation dom- 
ination occurred before the kinetic energy reached the 
potential energy in 0, the universe would return to an 
inflationary regime from which it would never recover. 
Our model has the same requirement. We will write a cn( j 
as the scale- factor when reheating ends (t ~ T^ 1 ) and 
the radiation redshifts as an unsourced relativistic fluid. 

The kinetic energy approaches the potential energy at 
a scalefactor a* , given by [l[ 

— < — m 4/3 Hm pi /M). (46) 
The energy in the <p field is 

p ^X m ^r x t m f(^) 6 (^)\ (47) 

\ a cnd/ V a / 

and the energy in the radiation is given in Eq. ([44)) . The 
ratio is of order unity when the kinetic energy reaches 
the energy in the radiation, at a scale- factor of a r . We 
require a r < a* to ensure radiation dominated expansion 
begins. Then, 

which is easily satisfied for the values of the parameters 
we have been considering, but does place a weak upper 
bound on A x for a given choice of g. 

Importantly, the energy in the potential of (f> varies 
only very slowly. The exact scale- factor a cn d < o>t < a*, 
where radiation transition occurs, matters very little in 
terms of its future evolution: the variation of the field's 
potential energy in the range in which this can happen is 
negligible, as the field is moving so slowly. The same evo- 
lution outlined in [l[ therefore takes place, with (f> mim- 
icking a cosmological constant and today's dark energy 
density obtained with M ~ 10 6 GeV . 

VII. DOMAIN WALLS 

Typically in models with a symmetry breaking term, 
domain walls inevitably form via the Kibble mechanism 
[l8| . Parts of the universe causally separated have no way 
of being correlated. When the x becomes unstable, 



in some regions the field will move to positive values, 
and in other regions to negative values. The domain walls 
created by the smooth transition between these values are 
generally a serious problem in many cosmological models, 
as they can come to dominant the energy density of the 
universe [l9j . 

Fortunately, and unlike in many models with such an 
occurrence, the symmetry is restored once <fi > m. Any 
walls formed will then dissolve. However, this could leave 
some effect upon the amount of reheating taking place 
within the regions of the domain wall, potentially in- 
fluencing large scale structure formation. As the wall 
thickness is significantly smaller than a horizon, we do 
not expect this to be a large effect, but warrants further 
investigation. 

VIII. SUMMARY AND FUTURE DISCUSSION 

By introducing an additional coupling to a second, self- 
interacting scalar field, we can restore the traditional re- 
heating mechanisms to the scenario proposed by Peebles 
and Vilenkin [![. The symmetry breaking term leaves 
the field temporarily unstable, which allows it to gain 
significant amounts of energy. Using a phenomenological 
approach to reheating, we have calculated the evolution 
of the relativistic particles produced by the decay of this 
field. We then described a variety of constraints on the 
model parameters, ensuring that the additional field docs 
not interfere with the behaviour of the original inflaton. 
The range in temperature produced is quite narrow if the 
new field does not strongly self- interact, and the symme- 
try breaking scale does not reach the Planck scale. But 
with a self coupling of the order of A x > 0.01, one has 

gl /4 T RH ~ 10 11 - 10 14 GeV. This allows us to suppress 
the contribution of the gravitational waves at the time of 
BBN well below the upper limit. 

With a model that tries to explain the evolution of the 
Universe from early inflation to today dark energy dom- 
ination, one of the issues to explain is the origin of the 
baryon asymmetry [20(. A possibility would be sponta- 
neous baryogenesis through a derivative coupling of the 
inflaton field to a matter current [2lJ. In our set-up, due 
to the larger coupling of x to the light degrees of free- 
dom (compared with the gravitational couplings), ther- 
malization occurs promptly after the start of reheating. 
This opens up the possibility for example of having lep- 
togenesis during or after reheating, by the decay of the 
lightest right handed neutrino. The field x being quite 
heavy during reheating could decay into the lightest right 
handed neutrino, or alternatively the neutrinos could be 
thermally produced. 

The model could potentially form domain walls, often 
causing problems in similar models. In this variation, we 
note they are transient and do not interfere with future 
evolution of the universe. The possibility exists their 
effect upon the reheating temperature will appear as an 
imprint on large scale structure formation. 
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The analyses of the reheating mechanism proposed in 
this letter has been done at tree-level. But having cou- 
pled the quintessence field to another scalar field, with 
a strength g > 10~ 6 , the question arises about the sta- 
bility of the quintessential potential to quantum correc- 
tions. Here we argue that indeed these corrections are 
under control due to the decoupling theorem, and that 
the scenario is not spoiled by quantum corrections, but 
leave the explicit calculation for a future work. 

Quantum corrections with the x field running in the 
loops can give rise to a potentially large correction to 
the quintessence potential, with the one-loop correction 
given by [13]: 

-™ = ^I>*>*(l»^-±), (-) 

*=<j>,x v p 7 

where /i is the renormalization scale and mi(4>) the field 
dependent masses, with m x — gcj). So for the x ne ld we 
have a very heavy state m x ~ gm p \ ^> gM that is excited 
in a universe with an energy density p -C AM 4 . Given 
that we do not have enough energy to excite such a heavy 
states, physically we can expect that they decouple from 
the spectrum [23| , and their contribution to the effective 
potential is highly suppressed. However, the 1-loop ef- 
fective potential as given in Eq. (|4*9")) is computed using 
a mass independent renormalization scheme that does 
not take into account threshold effects. The decoupling 
would appear naturally when using instead a mass de- 



pendent renormalization scheme [24]. To deal with this 
problem when working with the effective potential one 
can use instead the "improved" effective potential [25[ 
by replacing all parameters (masses, couplings and field 
vevs) by their renormalized values in both the tree level 
and one loop potential, and imposing the physical con- 
dition that the potential does not depend on the renor- 
malization scale, dV(p)/d In nl^—^, =0. By choosing the 
renormalization scale p* below any heavy mass threshold 
in the model, heavy states decouple and the dependence 
on fi is minimised. We are then left with the tree-level 
potential again with all the parameters evaluated at p* , 

U(<t>, X ) = V(&M0O,A(/**)) 

+^^_ TO V)) + ^X 4 . (50) 

All possible quantum corrections due to the heavy states 
are then encoded in the running of the mass parameters 
and couplings through the renormalization group equa- 
tions, and are therefore expected to be under control. In 
future work, we plan to do a detailed analysis of quantum 
corrections and decoupling in this model. 
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